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Abstract 

The eigenvectors of the Hamiltionians of the XYZ Gaudin model are 
constructed by means of the algebraic Bethe Ansatz. The construction 
is based on the quasi-classical limit of the corresponding results for the 
inhomogeneous higher spin eight vertex model. 
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1 Introduction 



The XYZ Gaudin model was introduced by M. Gaudin H, |, |] as a quasiclassical 
limit of XYZ spin-1/2 chain. Gaudin noticed also that the former model can 
be generalized to any values of constituing spins. Whereas the spectrum and 
eigenfunctions of the XXX and XXZ variants of Gaudin model can easily be 
found via Bcthe ansatz, the general case encounters the same problems as in case 
of the original XYZ model. For the spin-1/2 XYZ model a proper generalization 
of Bethe ansatz was found in Q], see also [||. However, when passing to the 
Gaudin model, an additional problem arises: whereas it is easy to perform 
quasiclassical limit in the Bethe equations determining the spectrum of XYZ 
model, the expressions for the eigenvectors appear to behave singularly. Solution 
to this problem is the main subject of the present paper. 

We obtain the arbitrary spin XYZ Gaudin model as a quasiclassical limit 
of the inhomogeneous higher spin generalization of XYZ model introduced in 
[H 0' Hi- To construct the eigenvectors of XYZ model we use the algebraic 
version of Bethe ansatz proposed in for spin-1/2 and generalized to higher 
spins in 0, |^. The paper is concluded with describing a way to circumvent 
the divergencies in the quasiclassical limit and to get finite expressions for Bethe 
eigenvectors of Gaudin model. 



2 Inhomogeneous XYZ model 

First let us briefly recall the inhomogeneous XYZ spin chain model. The model 
is parametrized by an anisotropy parameter rj besides the elliptic modulus r, 
shifts of the spectral parameter Zn, and a spin at each site In- Later we will take 
the limit — > and recover the generating function of the integrals of motion of 
the XYZ Gaudin model as the leading term of the transfer matrix of this model 
as in j^. 

The inhomogeneous XYZ model is a quantum spin chain model defined on 
the Hilbert space H = (2)^=1 Ki, where Vn is the spin in representation space 
Goo""^ of the Sklyanin algebra |lo|, jl|. See Appendix for notations. 

The transfer matrix of the model, t(u), is a linear operator on Ti. defined as 
the trace of the monodromy matrix, T{u), in the context of the quantum inverse 
scattering method 

T{u) = ^^n{u)) ■^^n{u-zn)...L2{u-Z2)Li{u-zi){2.\) 

i{u) = trc2{T{u))^AN{u) + DN{u), (2.2) 
where the L operators, Ln{u), whose elements act non-trivially only on the n-th 
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component, Vn, are defined by (cf. ( [A.2| )) 

^-(«) - (;:;:,' tS) 

p^^-iS") = 1®...® l^p^-CS"")® 1®...® 1. (2.4) 

Because of the fundamental commutation relation, 

Ruiui - U2)Toi{ui)To2{u2) = To2{u2)Toi{ui)Ri2{ui - U2), (2.5) 

the transfer matrices commute with each other: 

[i{ui),i{u2)]^0. (2.6) 

The quantum determinant ( |p^ ) is defined by 

A{u) tri2 P,-2Toi{u - r])To2{u + tj), (2.7) 

where is a projector onto the space of anti-symmetric tensors. This operator 
commute with all elements of the Sklyanin algebra and acts on the Hilbert space 
H as a scalar multiplication: 

A{u)\n = A_(u-?7)A+(u + ?7), (2.8) 

eii(u - Zn) 
n— 1 ^ ^ 

3 XYZ Gaudin model 

We define the XYZ Gaudin model as a quasi-classical limit of the inhomogeneous 
XYZ model defined above. Let us examine the asymptotic behaviour of the 
operators in the previous section when rj tends to 0. The L operator, the 
monodromy matrix, the transfer matrix, the quantum determinant and the R 
matrix are expanded as 



Explicitly 



Ln{u) = 


1 - 


V2f^C.a{u)+0{lf), 


(3.1) 


T{u) = 


1 - 


^2rjT{u)+r]^T^^\u) + 0{f]^), 


(3.2) 


i{u) - 


1 - 


Vrf trr(2)(w) + 0(7^3), 


(3.3) 


A(^) = 


1 - 


VTf{trT^^\u)+AtrT(uf) + 0{ff), 


(3.4) 


R{u) = 


1 - 


- 2rir{u) +0{ri^). 


(3.5) 
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N 





£„(u) = Wa(u)p'!^S'' ® Cr", T{u) = ^ Cn{u ~ Z„). 



a=l 
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(Tio(u) 

a{u) 


cn 

= — (" = 
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Oio 


9io(u) 
9ii{u)' 


o-oo(w) 


dn 

= —m = 

sn 






a{u) 




0ii{u) 


aoi{u) 


1 






a{u) 


sn(it) 







Here denotes the spin £ representations of the Lie algebra sl{2), and 
wi{u) - 

W2{u) - 

wsiu) - 

where 9ab — ^ab(O), 0[i = d/du{9ii{u))\u=o, <J is Weierstrafi' sigma function and 
5" are generators of the Lie algebra sl{2): 

Here {a,b,c) denotes a cyclic permutation of (1,2,3). The commutation rela- 
tions of the £ and the monodromy operator are 

[£.oi{u),Co2{v)] = [r{u - v),Coiiu) + £.02iv)], (3.6) 
[Toi{u),To2{v)] = [r{u-v),To,{u)+To2{v)]. (3.7) 



The classical r matrix r is defined by (3.5), or explicitly by 

1 ^ 

a=l 



The commutation relation (3.6) is nothing but the quasi-classical limit of the 
fundamental commutation relation ( |2.5| ). 

We define the XYZ Gaudin model by specifying its generating function of 
integrals of motion as f(u) — trT{u)'^. It is obvious from (3.3) and (3.4) that 
f{u) can be expressed as follows (cf. [H): 

A{u) - i{u) - 1 = Vf(u) + 0{Tf). (3.8) 

Therefore we can expect that eigenvectors of t{u) can be constructed as a leading 
term of the 77 expansion of eigenvectors of i{u) which are found in 0, 
Essentially this is true, but we must be careful in taking the limit as we will see 
in the next section. 

Operator t{u) is explicitly written down as follows: 

N N 
t{u) = ^ p(w - Z„)C(C + 1) + I] H.n({u ~~ Zn) + Hq, (3.9) 
n=l n—1 

where and p are Weierstrafi' zeta and p functions and 

3 N 

Hn = 2 ^ ^u;,(z„-z™)5^5;^, ^ = 0, (3.10) 

m^n a— 1 n—1 
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N 3 / ^z„-z,„+(^5/2 



n— 1 a— 1 \ m^n 
N 3 



/2 

+ E ^a{Zn - Z„,) (C(z„ - Z„ + ^) - C(^)) S^^S:^) (3.11) 

are integrals of motion. Here a is 1,3,2 for a = 1,2,3 respectively, wi = 1, 
L02 — T, LOy, — 1 + T, Ca = p{uJa/'2.). We Omit the symbol for the sake of 
simplicity. 



4 Algebraic Bethe Ansatz for XYZ spin chain 

The algebraic Bethe Ansatz in the context of the quantum inverse scattering 
method is applied to the inhomogeneous XYZ spin chain model in the following 
way. (cf. [^, Q, [^), Hereafter we assume that the total spin ^totai ~ (■! + • ■ •+^Ar 
is equal to an integer M. 

First we introduce the matrix of the gauge transformation: 



A twisted monodromy matrix is defined by means of My. 



Tx,\'{u;v) 



A\,\'{u;v) Bx,x>{u;v) 
Cx.x'{u;v) Dx,x'{u;v)^ 

Mx{u + v)-^T{u)Mx'{u + v), (4.2) 



where i; is a parameter. There is a vector ujxiv) called a local pseudovacuum in 
the spin £ representation space Vi = Oqo^: 

cjx{v)=cjx{v;y):= Ylf=i ^lo (^y + ^ - ^ + (2j - 3£ - 1)?]) x 

X ^10 (y - ^ + ^ " (2j -3£- 1)77 



and a global pseudovacuum in Ti: 

n^iv) := Wa+4(£i+...+£„)^(w + zn) 



®t^X+4(t.i + ---+e„)rj{v + Z„) (8) . . . ® UJx+4£i7j{v + Zi). (4.3) 



4 



As shown in j|] and we can construct eigenvectors of the transfer matrix 
t{u) of the form 

'f,{wi,...,WM) = ^e2™$A+2a^(ii;i,...,u;M), (4.4) 

aGZ 

$a(wi, . . . , Wm;w) Bx+2r,,X-2riiwi)Bx+4:r,,X+4ni'^2) ■ ■ ■ 

■■■ B\+2Mri.\-2Mr^{wM)^]i'^^^'''{v), (4.5) 

where an integer u and the parameters wi, . . . ,wm satisfy a system of equations 
(Bethe equations) for j — 1, . . . , M: 



e 



A_ (wj) k=ii=^j - ^"fe - 2?7) ■ 

The eigenvalue t{u) of the transfer matrix t{u): 

i{u)'^^{wi, . . . , wm) = t{u)'^^{wi, . . . , wm), (4.7) 

satisfies 

t{u)q{u) = A+(u)g(u - 2r]) + A_{u)q{u + 27]), (4.8) 

where 

M 

g(u) =e— " J] ^nlw-w'n). (4.9) 

m— 1 

Note that the Bethe vector '^/^{wi, . . . , wm) is expressed as a Fourier series, 
and the convergence of such series is not a priori known. In fact when rj = r' jr 
is a rational number, this Fourier series diverges and we must replace it by a 
sum with r terms. Hence we cannot take the limit 77 ^ of the above results 
naively to obtain eigenvectors of the XYZ Gaudin model. 



5 Algebraic Bethe Ansatz for XYZ Gaudin model 

In order to obt ain an eigenvector of the XYZ Gaudin model, we bypass the 
divergent series (4.4) in the following way. At an intermediate stage of derivation 
of the Bethe equations, we have the formula: 

i(u)^\(w;v) = A(u;w;77)$A-27,(w;w) + A(u; zZ;; -r/)$A+2>,(u'; «) + 

M M 

+ J2^jiu;w;Tj)^x-2r, {wj (u) ; v) + ^ A^^ (m; uJ; ~?7)$a+2,, (wj {u);vl^.l) 
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where we abbreviated (wi, . . . ,wm) to w and {wi, . . . ,Wj-i,u,Wj+i, 
to 'Wj{u), and A and are defined by: 



M 



^\ Oll{u - Wm) 
m— 1 ^ ^ 



A^(?x;w;;ry) = 0ii(-2r,). 



n 



■ll{Wj - Wk) 



■ ,Wm) 

(5.2) 
(5.3) 
(5.4) 



Since all of these formulas are holomorphic around rj = 0, we can expand them 
as Taylor series in 77. For example, as M\{u) docs not depend on 77 and the mon- 
odromy matrix T{u) has the expansion (3.2), the twisted monodromy matrix 

T\+2mr,,\-2m7ji'u;v) has the expansion. 



{u;v) 



1 
1 



oiv). 



(5.5) 



Therefore (l,2)-element of the twisted monodromy matrix has an expansion as 
follows: 

Bx+2mi-,,X-2mr){u\ v) = rjBx^miu, v) + 0(7?^). (5.6) 

The pseudovacuum vector defined by (^^) is expanded as 



J7^-2^^''(.)=17^^-^^(°)(.)+,;17^^-^'W(.) + 



Combining (5.6) and (5.7), we obtain 

$a(w;w) = ?7*'^((/)a(w;w) + 0(77)). 



(5.7) 



(5.8) 



It is possible to write down the complicated definition of (t>x{w; v) explicitly, but 
it is not necessary to our purpose here. Expanding A (5.2) and Aj (5.4) as 

A{u]w-r]) = l + r]A'^^\u]w)+?fA^^'>{u;w) + 0{Tj^), (5.9) 
A^{u-w-r]) = T]A^''^\u;w)+r]^A^''^\u;w) + 0{7j^), (5.10) 

we obtain the following equation for the action of leading terms of t{u) on 
(j)x{u;w;v) from (O) and (p.3[): 



trT'^^\u)(f)xiw;v) 

- - 4A(i) (u; w)-^ + 2A(2) {u; w)^ (w; v) 

1=1 \ / 
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Instead of the Fourier series (4.4), we take a Fourier transformation of 
b\{w; v) as a candidate of an eigenvector of trT(^)(u) (and hence of t{u)): 

»i 



Mw) - / e^'''^cj)x{w;v)dX. (5.12) 
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Here v is an integer to be determined. Since <j)\{w]v) is periodic with respect 
to A with period 2, operator trT*^^) acts on tp^ as follows by virtue of ( ^.11 ): 

trr(2)(u)V'i.(w) = {A{mvY + AmvK^'^\u]w) +2k'^'^\u]w))i}y{w) 
M 

+ (unwanted) j, (5.13) 

Here "unwanted terms" are defined by 

(unwanted)j = j e^'"^ (^^AKf^\u\ + 2kf'^\u- vo)^ (t,x{wj{u)-v) dX 



e 



1 I ^1 &lliWj~Zn) SlliWj - Wk) 

X 4Af^\u;w)(l)x{wj{u);v)dX, (5.14) 



where ' denotes the derivative. Thus, if (...) in the last expression of (5.14) 
vanishes, namely, if 

holds for any j = 1,...,M, then ( 5.13| ) implies that ipi,{wi, . . . ,wm) is an 
eigenvector of tr T^^^ . 

The corresponding eigenvalue of f{u) can be computed from the expan- 
sion {p.9\j, but a simpler way is the following: Observe that the eigenvalue of 
trT(2)(^) for iP^iw) 

A{mvf + AmvK^^\u\ w) + 2A(^)(m; w) 



is equal to the coefhcient at rf of the right hand side of (4.8) in which we replace 
{wj} by a solution to ( |3.15 ) this time. Thus we can derive the formula for an 
eigenvalue of f{u) corresponding to ipy by extracting terms of order rf in (2 
and Q (cf. @): 



q"{u) - 2Z{u)q'{u) + {Z{uf - Z{u)')q{u) = T{u)q{u), (5.16) 

or 

t{u) = ixiu) - Z{u)f + ^ixiu) - Z{u)), (5.17) 

du 
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where 

This result is an eUiptic analogue of the spectrum of the XXX Gaudin model 

i, i, i, §• 

6 Comments 

We applied the algebraic Bethe Ansatz to the XYZ Gaudin model and con- 
structed eigenvectors of the generating function of integrals of motion, f(w). 
Several questions need, however, further investigation. 

First, it would be preferable to have a closed treatment of the Gaudin model, 
completely independent of the original XYZ model. The difhculty is that the 
expression for (j)\{'W] v) obtained by differentiation of <^\{w; v) is rather volumi- 
nous. 

Another weak point is that we must assume that the total spin ^totai must 
be an integer. This obstruction comes from the same assumption in the Bethe 
Ansatz of the XYZ type spin chains [|| , |^ . 

It is also an open question if we can obtain all the eigenvectors whithin our 
approach. However, this is quite common situation in the application of the 
algebraic Bethe Ansatz known as completeness problem. 

We expect that these problems would be overcome within an alternative 
approach to the model known as Separation of Variables, The functional 



equation (5.16) should be interpreted then as a separated equation. 



A Review of the Sklyanin algebra 

In this appendix we recall several facts on the Sklyanin algebra and its repre- 
sentations from and We use the notation from jl^ for theta functions: 

Oab{z] t) = X! ''^P {'^^ (i ") ^ ^'^^ (i ^ ") ( ^ ^ ^ 

riGZ ^ ^ 

where r is a complex number in the upper half plain. 

The Sklyanin algebra, Ut.ii{sI{2)) is generated by four generators 5*^, 5*^, 5*^, 
5*^, satisfying the following relations: 

i?12(w - u)Loi(w)Lo2(w) = L()2{v)Lqi{u)Ri2(u - v). (A.l) 

Here u, v are complex parameters, the L operator, L{u), is defined by 

3 

L{u) ^ y (A.2) 

a=0 
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2^?n(77;T)' ^ ^ ^ 20n(u; t)0io(??; r) ' 



R{u) is Baxter's R matrix defined by 

3 

R(u) = Wa {u)a'' ® a^ l^i^(u) W^{u + + ry). (A.3) 

and indices {0, 1,2} denote the spaces on which operators act non-trivially: for 
example, 

3 3 
i?i2(w) Loi{u)^J2^aiu)S'' f^a" (8)1. 

a=0 a=a 

The above relation ( [A.l| ) contains u and f as parameters, but the commu- 
tation relations among S"" (a = 0, . . . , 3) do not depend on them: 

[5", = -iJa.AS", [S'',S''U = i[S", S-'U, (A.4) 

where (a, (3, 7) stands for any cyclic permutation of (1,2,3), [A. B]± = ABztBA, 
and J„,0 = (W^ - Wl)/{W^ - Wi), i.e., 

"^12 — 7i — 7 ^97^ — 7 TT' "^23 — 7; — 7 TTT — 7 r:T, ^31 — ^ - 



9oo(^;T)20io(r7;T)2' 6oa{mrY0^i{mrY' 9oiimrrdioimry 

The spin £ representation of the Sklyanin algebra, 
/:!7,,^(sZ(2))^Endc(e«) 
is defined as follows: The representation space is 

©00+ = {/(^) I /(^ + 1) = /(-^) = /(^), /(^ + ^) = CXp-«..(2.+r) ^(^)i^ 

It is easy to see that dim0QQ+ — 21+1. The generators of the algebra act on 
this space as difference operators: 

9ii[2z;t) 

where 

soiz) = 6*11(77; T)6'ii(2z;r), si{z) = 6*10(77; t)6'io(2z; t), 

S2iz) = i6'oo(?7;T)6'oo(22:;r), 53(2) = 6*01(77; t)6'oi(2z; r). 

These representations reduce to the usual spin £ representations of U{sl{2)) for 
(r/^0). 
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